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ABSTRACT 

In recent work, we have shown that it is possible to hnk quantitatively many aspects of damped 
Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of 
galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we 
consider the expected properties of intrinsic Lyman alpha absorbers seen in the spectra of 
high redshift {z > 2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with 
the death of massive stars, their afterglows provide insights into otherwise unprobed regions 
of protogalactic objects, but detailed physical interpretations are currently embryonic. 

We find that median impact parameters (measured from the potential minimum) are 
approximately 1 kpc for GRBs compared with 4kpc for QSO-DLAs. However, an equally 
important difference is that GRB-DLAs are predominantly associated with halos of mass 
10^" < M^iJMq < 10^^, an order of magnitude larger than the hosts of QSO-DLAs. Ac- 
cordingly, there are differences in the stellar properties of hosts. For instance mean star forma- 
tion rates are higher: (A/^) ~ 10 Mgyr^i for GRB-DLAs compared with (M^,) ~ IMgyr'^ 
for QSO-DLAs. 

Our simulations accurately predict the form of the GRB-DLA H I column density dis- 
tribution, producing quantitative agreement for A'hi > 10^^ cm^^, but they somewhat un- 
derpredict the incidence of low column densities A^hi < 10^^ cm^^. This is reflected in our 
estimate of the ionizing photon escape fraction, /osc — 1%, which is lower than the obser- 
vational GRB-derived escape fraction (2%). Line-of-sight neutral gas metallicities predicted 
by our simulations (10^^ < Z/Zq < 1) are consistent with the modest observational con- 
straints. Because of large internal dispersions in gas metallicities, this agreement is not sig- 
nificantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs 
(Zi, < Zq/3), confounding claims that the observed metallicity of GRB-DLAs poses a severe 
challenge to current GRB models. 
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1 INTRODUCTION 

Observations of transient Gamma Ray Burst (GRB) afterglows 
have recently been established as a promising new window on the 
high redshift Universe. In particular the coupling of rapid response 
observing opportunities on large telescope facilities with the de- 
tection alerts provided by the Swift satellite ( Gehrels et al. 2004 1 
have produced a modest but growing sample of high quality data. 
The GRBs are assumed (in the case of long duration GRBs which 



* Email: apontzen@ast.cam.ac.uk 



have brighter optical afterglows and will be the focus of the present 
work) to be associated with the death of Wolf-Rayet stars ( |Mac-| 
[Fadyen & Woosley|1999| ). 

Intrinsic absorption features in the GRB afterglow spectra are 
of particular interest because they probe the gaseous environment 
of the progenitor. The physical state of the gas may be quite dif- 
ferent from that of typical (cross-section weighted) intervening ab- 
sorbers. Dense neutral hydrogen giving rise to intervening absorp- 
tion with column densities TVhi exceeding 2 x 10^" cm^^ is known 
as a damped Lyman alpha system (DLA); the Hi absorption and 
its associated metal line features are commonly studied in the spec- 
tra of quasars (e.g. [Wolfe et al.||2005] l. The focus of the present 
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|et al.|2006); whereas the median QSO-DLA metalhcity is around 
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Figure 1. A sketch of relevant differences between typical QSO (dashed 
vertical line) and GRB afterglow (solid vertical line) sightlines through a 
protogalaxy at high redshifts (z ^ 3), after [Prochaska et al.H2007) . As- 
suming GRBs are associated with the death of massive young stars, this 
class of object arises inside dense H II regions (white circles) within the 
the most rapidly star forming regions of the neutral medium (solid black 
regions); whereas QSO sightlines are more likely to intersect the neutral 
medium at a larger impact parameter In our simulations, single protogalac- 
tic halos (particularly those with A/vir > 10^"'^ Mq) can host multiple 
neutral star-forming regions; thus GRB sightlines may also intersect a sepa- 
rate neutral zone which has little physical association with their origin. The 
local stellar UV, presumed largely absorbed within the H II regions by our 
model (Section [2. 1) , has an effect on GRB-DLAs but appears unimportant 
in understanding the QSO-DLAs. However the incoming background inter- 
galactic mean field has an impact on both classes of object (Section 
Typical temperatures and densities for the regions of interest are indicated 
(based on observations of our own Galaxy). 



paper will be the analogue of DLAs seen in the intrinsic spectra 
of GRB afterglows: neutral hydrogen absorption at the emission 
redshift of the GRB itself, by analogy termed a "GRB-DLA'Q Be- 
cause GRB-DLAs are physically associated with the GRB hosts, 
their nature may be quite different from intervening DLAs studied 
in quasar spectra (QSO-DLAs); part of the aim of this paper is to 
recap known observational differences and link them to the under- 
lying protogalactic populations. 

Indeed, a statistical comparison of QSO-DLAs with GRB- 
DLAs immediately confirms that the populations are distinct. 
While the frequency of quasar absorption systems decays monoton- 
ically with increasing neutral hydrogen column density (e.g. |Wolfel 
|et aI.|2005[[Prochaska et al.|2005[ l, the observed frequency of GRB- 
DLAs grows to a peak at A'^hi — 10^^ cm^^ (a column density 
exceeding the strongest QSO-DLA ever observed; Jakobsson et al. 
2006). Metallicities of the H I-associated gas, although challeng- 
ing to measure in the transient GRB-DLA afterglows, are generally 
in the range 1/100 < Z/Zq < 1, with a median somewhat ex- 
ceeding Zq/10 (e.g. Savaglio|2006 Prochaska et al.|2007 Fynbo 



The physical interpretation of the GRB-DLA population is 
maturing, but is not firmly established. It seems clear that, quali- 
tatively, an association of GRB events with the death of massive, 
short-lived stars leads to the correct expectation of higher column 
densities (to ensure high star formation rates) and higher metallic- 
ities. This has led to the schematic picture (e.g. [Prochaska et al.| 
2007 ) that GRB-DLAs preferentially probe the innermost regions 
of rapidly star forming protogalactic halos (Figure[TJ, while QSO- 
DLAs intersect the outer regions of such systems. The neutral gas 
probed by the GRB-DLAs is not thought to be in the immediate 
circumstellar environment but rather at distances exceeding at least 
10 parsecs. This follows because, under realistic assumptions, the 
UV emission from the afterglow itself is capable of ionising any 
local neutral material within just an hour of the photons from the 
initial burst arriving (Proch aska et al.|2006j >. The non-local picture 
is reaffirmed by the non-detection of time variations in the neutral 
columns ( |Pema & Loeb|1998[[Mirabal et al.|2002[ l. In addition, ap- 
proximate upper limits on distances of absorbing gas can be derived 
from UV pumping of fine-structure transitions, which suggest the 
responsible neutral gas is at a distance of at least 1 00 pc ([Prochaska 
et al.''2006) and perhaps as much as ~ 1 kpc ( [Vreeswijk et al. 



For simplicity this term is applied to all intrinsic GRB absorption systems 



even though 20% would not exceed the traditional DLA threshold. 



2007; D'E lia et al.[|2009) . Finally, dust depletion patterns seem 
more consistent with those seen in the warm neutral ISM of our 
Galaxy, rather than in cold dense molecular clouds jSavaglio|2006[ l. 

It is possible that, in addition to the GRB-DLAs probing par- 
ticular regions of halos, the population of objects giving rise to 
GRB-DLAs are on average more massive than those giving rise to 
QSO-DLAs. This will arise if the cold gas cross-section (relevant 
for QSO-DLAs) has a shallower than linear relationship with the 
UV luminosity (assumed relevant for GRB-DLAs); such a picture 
has been successfully employed by |Fynbo et al.[ ( [2008| l to repro- 
duce differences in metallicities of the absorber populations. 

As well as providing insights into individual objects, GRB- 
DLAs offer an unusually direct way to probe the escape of ionizing 
photons from the star forming regions with which they are associ- 
ated. On the assumption that the GRB rate traces the star formation 
accurately, one may construct a global escape fraction (/esc) — 2% 
by modest extrapolation from the GRB-DLA column density dis- 
tribution (Chen et al.|2007[ Fynbo et al.|2009) . Conversely recent 
simulations by |Gnedin et al.[ ( [2008| l, aimed at constraining escape 
fraction statistics, have incidentally been able to account for the 
GRB-DLA distribution (albeit with some mismatch at high column 
densities). 

The emphasis in our work is reversed: we focus on statistics 
of GRB-DLAs between 2 < z < 4, 90% of which have column 
densities exceeding 10^^ cm~^. We use high resolution hydrody- 
namical simulations with physically motivated star formation feed- 
back algorithms ( [Stinson e t al. 2006) to generate a catalogue of 
GRB-DLA sightlines. This catalogue can be compared to current 
observational constraints, while the association of each simulated 
sightline with a specific protogalactic object allows us to predict 
the relationship between GRB-DLAs and other known classes of 
objects (see also [Fynbo et al.[[2008^ . We can regard such associ- 
ations with some confidence since, after calibrating two free pa- 
rameters in isolated disk galaxy tests, the simulations are known 
to reproduce realistic z — disk galaxies ( [Governato et al.|2007| 
[2009) , have realistic stellar mass-metallicity relations at low and 
high redshifts ( [Brooks et al.[ 2007; see also section 7.5 of |Maiolinol 
[et al.| 2008), and have produced a population of QSO-DLAs which 
are as realistic as current techniques permit (for a detailed compar- 
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ison to observations see IPontzen et al.|2008| henceforth P08). The 
treatment of local ionising sources we adopt is simple compared to 
detailed modelling employed in some ionising photon escape frac- 
tion studies (e.g. Razoum ov & Sommer-Larsen|2006[|Gnedin et al.| 
12008, Wise & Cen 2009 1. As a result, the statistics of systems with 

I rn ^ ' 

A'hi < 10 cm seem to be sensitive to assumptions in our ap- 
proximations, and corresponding caution will be exercised in our 
interpretation - but results appear robust for larger column densi- 
ties. 

The simulations employed and techniques applied in the 
present work are closely related to those of P08, and are discussed 
further in Section|2] By tracing sightlines from star forming regions 
we produce a set of simulated GRB-DLA statistics which we com- 
pare, where possible, with existing observations (Section[3]l. While 
both observational statistics and theoretical models are far less cer- 
tain for GRB-DLAs than QSO-DLAs, we are able to draw tentative 
conclusions from our present study in Section]?] 



2 SIMULATIONS AND TECHNIQUES 

In this Section, we describe the simulations employed and discuss 
the techniques required to build a cosmological sample of GRB- 
DLAs. 

The simulations considered here are similar to those employed 
by P08. In particular, they take advantage of the 'volume renor- 
malization' technique, in which the dynamic range is expanded by 
computing a subregion of the full volume at high resolution l |Katz| 
|& White|I99"3) : all results are taken from the high resolution zone. 
Time evolution is computed with the SPH code GASOLINE ( |Wads-| 
[ley et al.|2 004 ) which implements gas cooling, the effects of a cos- 
mological UV background ( [Haardt & Madau|I996^ with a local ap- 
proximate self-shielding algorithm (see section 5.4 of P08), and the 
effects of star formation based on the description of jStinson et al.| 
( |2006) . Supernova ejecta abundances for oxygen and iron are cal- 
culated in the simulation following the yields of [Thielemann et aT] 
( |1986^ and |Weaver & Woosley| ( [T993] l; after being injected into the 
interstellar gas these metals are advected with the SPH particles, 
although we now incorporate a correction for sub-resolution tur- 
bulent mixing which can transfer metals between Lagrangian fluid 
elements (see |Shen et al.|2009j l. Gravitational softening is constant 
in physical units for z < 8 and evolves comovingly for z > 8; its 
final fixed value is given for each simulation in Table [T] For a full 
description see section 2 of P08. 

The particular regions in use for the present work are listed 
in Table[T| These mirror closely the regions used by P08, although 
none are in fact identical. The first three are selected to form (by 
2 = 0) individual objects and their immediate environments: 'DG' 
is a high resolution run which, at 2: = 0, forms a dwarf galaxy 
(A'/vir ^ 3 X 10^° Mq) at 2 = 0; 'MG' forms a Milky Way-like 
galaxy (Mvir ~ 7 x lO" Mq) and 'CL forms a cluster (Mvir ^ 
1.1 X 10"Mq). 'C50' is a full cosmological box of side 50 Mpc 
(comoving). From each box we select, at 2 = 3, the set of halos 
having greater than 3 500 dark matter particles - an empirical limit 
which we established ensured convergence of our results including 
stellar populations (see also P08; Brooks et al.||2007|f . The 'CL 
simulation, combining high resolution with large volume, contains 
86 halos satisfying this stringent criterion whereas the three other 
simulations have only ~ 15 each. For this reason, we selected a 
random subsample of 15 halos from the 'CL simulation for our 
statistical results which are therefore derived from 60 halos of high 
resolution, evenly distributed with 10** < Mvir/M© < lO"'^. We 



Tag 


max Afp,gas 


Ai'p,DM 


e/kpc 


Usable Vol (comoving) 


DG 




1O*-6M0 


0.12 


4 Mpc^ 


MG 


1Q5.6 


106.0 


0.34 


70 Mpc^ 


CL 


1Q5.6 


106.0 


0.34 


600 Mpc^ 


C50 


107.9 


108.4 


2.00 


125 000 Mpc^ 



Table 1. The simulations used in this work. The first column is the tag 
which we use to refer to each simulation. For all except "C50", a subsam- 
ple of the full box is simulated in high resolution; no results are taken from 
outside this region. The second and thii'd columns refer respectively to the 
maximum gas and dark matter particle masses within the region, the fourth 
to the gravitational softening length (in physical units) and the final column 
gives the comoving volume of the region. The separate boxes are generated 
from entirely different sets of initial conditions; the "DG" and "MG" sim- 
ulations are designed to form respectively a dwarf and approximately L^, 
galaxy at z = while the "CL" region forms a large cfuster of galaxies at 

2 = 0. 

discuss the impacts of the finite mass range of resolved halos in 
Sectionim 

As always with 'zoomed' simulations, one must be as care- 
ful as possible to test whether the arbitrary choice of region may 
have affected results through environmental dependences; we take 
agreement of statistics between boxes (where their resolved halo 
mass ranges overlap) as evidence that no severe problems are to be 
expected. In particular our present simulations include at 2 = 3 
an over-dense proto-cluster region ('CL) which resolves the same 
range of halo masses as the more typical 'MG' box. There is no 
evidence for meaningful disagreement between the properties of 
halos within these boxes (see also |Crain et al.|2009[ l. However, as 
in P08, we must accept in principle the possibility that we miss an 
environmental effect which has an impact on statistics - this is to 
be regarded as a trade-off against resolving such a large range of 
halo masses with relatively uniform resolution. 

While the regions are from a different realization of the initial 
power spectrum than those presented in P08, the physics remains 
almost identical except for two notable differences. Firstly we have 
used the self-shielding approximation, described in section 5.4 of 
P08, in all our simulations. This has the effect of slightly raising 
cold gas fractions and lowering dwarf galaxy metallicities (the lat- 
ter because the lower pressure of the ambient ISM makes our feed- 
back algorithm somewhat more effective); see P08 for details. Sec- 
ondly, we have implemented a turbulent metal diffusion algorithm 
as described at the start of this section. 

Because of the long lead-time involved in running the simu- 
lations considered, the cosmological parameters employed (Qm = 
0.24, = 0.76, D.t = 0.042, h = Ho/{WO kms'^) = 0.73, 
(jg = 0.77, Us — 0.95) are close to, but not fully, consistent 
with the latest constraints from large scale observational results 
(e.g. [Dunkley et aLl|2009^ . As we have previously argued, we do 
not believe small differences in the adopted cosmology will impact 
strongly on the results (see section 5.3 of P08). 

2.1 Approximate Radiative Transfer 

Understanding the effect of radiation on the ionization state of the 
gas is crucial to obtaining realistic absorption statistics. The sim- 
ulations incorporate a diffuse cosmological UV background based 
on |Haardt & Madaul ( |1996l l, uniform except in dense regions for 
which a local self-shielding algorithm, described in Section 5.4 of 
P08, is applied. We post-process the output with a simple equilib- 



© 2009 RAS, MNRAS 000, 000-000 



4 A. Pontzen et al. 



rium radiative transfer algoritiim, described in Section 3.1 of P08, 
to improve our estimate of the self-siiielding effect. We previously 
argued that the correction to this picture from local sources, which 
were not included, would be small. 

However, in GRB-DLAs, where the sightline necessarily 
passes through highly ionized bubbles around hot stars, local ra- 
diation sources are likely to have a stronger impact on results 
(these hot bubbles presumably having negligible cross-section to a 
QSO sightline). Observationally, fine-structure transitions in GRB- 
DLA absorbers indeed suggest the presence of a local UV field 
( |Prochaska et al.|2007) . 

A study of the local ISM shows that incorporating the stel- 
lar ionizing sources into a radiative transfer scheme using the raw 
density field from the simulations will overestimate the effect of 
local sources. This is because young stars are initially surrounded 
by a dense molecular cloud, which is rapidly photoevaporated to 
form an Hll region. Such regions have high densities, typically 
between 10 and 100 cm^'' but sometimes as high as 10* cm^"*. 
Hence their recombination time is short and they are able to absorb 
the vast majority of the UV radiation produced within them. This 
micro-physics occurs on parsec scales, below the resolution limit 
of our simulations. 

Assuming typical temperatures and densities for H II regions 
(T = 10* K, riH = 30cm~'^), one may estimate the hydrogen 
mass required to fully absorb all UV photons produced, using an 
on-the-spot approximation to account for photons emitted by re- 
combination to the ground state. Using STARBtJRST99 ( [Vazquez] 
|& Leitherer|2005l l emission models we find 



Mhi^hii 1.5 X IO^Mq 



IMoyr- 



30 cm 



(1) 



of neutral hydrogen would be ionized by the local sources, where 
Mi, is the star formation rate, the final factor is fixed by assumption 
in our calculations and a |Kroupa et al.| ( |T993^ initial mass function 
is assumed for consistency with the simulations. 

To assess the impact of removing the corresponding neutral 
material, we calculate the total star formation rate per halo and ex- 
pand ionized spheres around all young stellar particles (< lO^yr) 
until the total mass converted from H I to H II is given by l[T](. Be- 
cause the star forming threshold density (0.1 cm~^) in the simula- 
tions are lower than true H II region densities, the volume affected 
by this algorithm is rather larger than the physical sizes of H II re- 
gions would suggest, but by construction the mass affected is accu- 
rate. This should therefore reduce line-of-sight column densities by 
an appropriate fractior]^ We verified that the impact on the quasar 
DLA statistics is, as expected, insignificant, but the effect on GRB- 
DLAs is non-negligible (Section|3j. 

Our approach ignores the effect of dust which is capable of 
absorbing a significant fraction of ionizing radiation ( [Osterbrock &| 
|Ferland|2006^ ; since the quantity of dust in high redshift Hll re- 
gions is hard to assess, we merely admit that the magnitude of our 
correction could be reduced by including dust (see also Section|4]l. 
We further fail to account for molecular clouds, which could be 
important in understanding DLAs and GRB-DLAs (e.g. |Schaye| 
|2001| l. The small fraction of detectable H2 in most QSO-DLAs 
(|Ledoux et al.||2003,l and lack of H2 in many GRB-DLAs ([Tum 



|linson et ah 2007; although see |Fynbo et al.| 2006; |Prochaska et al 



^ There is some evidence to suggest that the GRBs themselves ionise mate- 
rial along their line-of-sight (see Introduction); with our current knowledge 
this is hard to model in detail, but we discuss possible consequences for our 
results in Section|4|. 



2009) does not imply H2 is unimportant, as large masses could 
conceivably reside in small clouds with negligible cross-section. 
|Ledoux et al.| l |2009^ also suggested that non-detections in most 
GRB spectra are expected given the low metallicities of typical sys- 
tems observed with high resolution spectroscopy. Overall the frac- 
tion of gas in H2 clouds is, at present, not calculable from our simu- 
lations and so we are forced to proceed on the assumption that later 
inclusion will make only a small correction to our results. Further 
discussion of the inherent problems in modelling the small-scale 
ISM and its impact on our results is given in Section |4] 



2.2 Sightline Generation 

The theoretical and observational picture of the precise origins of 
long GRBs is complex (e.g. [Woosley & Bloom||2006 l and refer- 
ences therein). We will suppose the events to be associated with the 
death of massive stars (~ 25M0); such stars have a main sequence 
lifetime of just a few million years, far smaller than the dynami- 
cal time of host halos (e.g. for the adopted cosmology at 2; = 3, 
(200 GpcY^^'^ ~ 7 X 10* yr). For the purposes of this work, long 
GRBs are therefore assumed to arise in a given region with prob- 
ability proportional to the star formation rate of that region. We 
choose the star particles younger than 5 x lO^yr in each halo as 
launch sites for our tracer sightlines in order to reflect this assumed 
association. Clearly the maximum age exceeds by an order of mag- 
nitude the actual lifetime of the GRB progenitors, but it is required 
to guarantee at least 200 candidate launch sites in each star forming 
halo (in tum required for reliability of our resampling mechanism 
described in Section [23) l. We verified that reducing cut-off times 
to lO^yr (the minimum possible at our resolutions) did not sig- 
nificantly affect our results with the exception of the low column 

- see Section '. 



3.2 



density distribution (< 10 cm~ ) - 

Additionally, we considered the effect of a metallicity ceiling 
for progenitor stars; such a limit is motivated by GRB coUapsar 
models which require mass loss over the progenitor Wolf Rayet 
stellar lifetime to be small ( [MacFadyen & Woosley |[r999[ l. How- 
ever, we found its effect - even on the measured column metallici- 
ties - to be minor. A discussion is given in Section[T6] 

The sightline properties are generated using the code de- 
scribed in P08, performing a full SPH line-of-sight integration to 
determine the column densities of neutral hydrogen and neutral- 
phase metals, i.e. the sightline metallicity is defined to be 



ghtlino 



dxZ(x) nm{x) 



(2) 



where x is the distance from the GRB along the line of sight and 
Z(x) is the metal mass fraction at that point. Equation |2| makes 
the implicit assumption that observers measure metal ions which 
are perfectly couplecj^to the H I (e.g., for silicon, that nsi 11 / nsi — 
"-Hi/nn). 

We use the oxygen mass abundance reported by the simu- 
lations as our metallicity tracer, using the reference value 12 + 
logj^Q(no/nH)0 = 8.69 (from [Lodders|2003| to normalize to a 
solar metallicity scale. The choice of oxygen reflects its status as 
the most abundant metal by mass; differences between elements 



This assumption was also made for generating Sill profiles in P08; 
however contrary to the statement in [Tescari et al.| 2009, we did not as- 
sume a fixed metal abundance when generating Si II profiles. We also used 
Cloudy models to show that the correction to the perfect coupling model 
is small (P08 section 5.2). 
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Figure 2. The product of the star formation rate and halo mass function 
for each qualifying halo (see Section |2] for resolution criteria) in our DG 
(crosses), MG (circles), CL (tripods) and C50 (plus symbols) runs. This 
quantity is proportional to the contribution of a given halo mass to the over- 
all cosmological ensemble of GRBs and therefore GRB-DLAs. A fit to the 
equivalent quantity for quasar DLAs is shown as a dashed line, based on the 
locus of points in figure 5 of P()8. Halos with no star formation are grouped 
together at the bottom of the figure. 



are small relative to the overall spread of metallicities under con- 
sideration. 



2.3 Resampling 

The set of halos extracted from all boxes is used to generate a cos- 
mological sample in a directly analogous fashion to that described 
by P08. However, instead of using a gaseous cross-section we ar- 
range that the probability of selecting a sightline from a particu- 
lar halo of virial mass Mh is proportional to the product of the 
halo mass function f{Mh) and the star formation rate of that halo 
Mi,ji. We bin our ensemble of halos in virial mass; the probability 
of a given GRB originating in the halo h from mass bin i{h) is then 
proportional to Wh, 



Wh = 



(3) 



where Fi is the halo mass function integrated over the mass range 
of bin i (P08 equation 4) and rii is the total number of halos from 
our simulations in that bin. All equations from P08 follow identi- 
cally with our replacement definition ([3]l for the halo weight Wh- 



3 RESULTS 

3.1 Origins of GRB-DLAs 

Because GRB-DLAs are intrinsic to the galaxy hosting the GRB 
event itself, the cosmological regions probed by such absorbers are 
not necessarily the same as those probed by the population of in- 
tervening QSO-DLAs. By assumption, our method triggers GRBs 
in a given region proportionally to that region's contribution to the 
cosmological star formation rate. This can be split into two factors: 



the relative contribution of halos of different masses, and the posi- 
tion of the star forming regions within each halo. In this section we 
will consider each of these two contributing effects in turn. 

Figure [2] shows the product of the halo mass and star forma- 
tion rate (/(M^) Af*.h) for individual halos in each of our boxes; 
as described in Section [ZJ] this relationship defines the relative 
contribution of halos of the specified virial mass. We find that a 
broad range of halos contribute to the overall population of GRB 
hosts, peaking at ~ IO^^Mq. The distribution is in good agree- 
ment between our four separate boxes in their regions of overlap. 
In particular the proto-cluster (CL) environment shows no system- 
atic differences from the disk galaxy progenitor region (MG). 

For masses Afvir < lO^'^M©, the contribution tails off as star 
formation rates are stifled by the stellar feedback and UV field (e.g. 
[Read et al.|2006| P08). The steep drop at low masses in Figure |2] 
shows that our low-mass cutoff of Afvir ~ 10* M0 should not lead 
to any systematic effects. For Afvir > 10^^'^ Mq the GRB rate 
contribution begins to be suppressed by the downturn in the cos- 
mological halo mass function; however, somewhat unsatisfactorily, 
our most massive halo has A/vir = 10^'^'^Mq which means that 
the volume of our simulations is insufficient to trace objects until 
they are so rare as to be cosmologically insignificant. We checked 
for systematic effects by artificially cutting off our distribution at 
IO^^Mq, finding that no results were severely effected, but a satis- 
factory resolution would require considerably larger simulations to 
produce a full population of rare A/vir > 10^^'^M© objects. 

We have overplotted (dotted line in Figure[2j a fit to the equiv- 
alent halo-weighting quantity for QSO-DLAs, o-dla/(A/) where 
o"DLA is the cross-section of the halo for column densities of neutral 
hydrogen A^hi > 10^" '' cm~^ from P08. The dominant contribu- 
tion to the QSO-DLAs can be seen to arise in halos approximately 
ten times less massive than for GRB-DLAs; the mechanisms de- 
termining the upper and lower limits are similar in both cases, but 
because ctdla scales less steeply than Af* (ctdla ~ Afvir, P08 fig- 
ure 4; Af* ~ A4if,P08 figure 13), the latter distribution is skewed 
towards higher masses. We will comment in section [33] that this 
has important consequences for the metallicity distribution of our 
simulated GRB-DLAs. 

We now turn to the location of GRB-DLAs within individual 
halos. Figure [3] shows the impact parameter h (projected distance 
of the DLA sightline to the minimum gravitational potential of the 
host halo) for GRB-DLAs (thick solid line) and QSO-DLAs (thick 
dashed line) in our current set of simulations. The former arise con- 
siderably closer to the halo centre (typically at ~ 1 kpc) than the 
latter (at ~ 4 kpc; although there is a spread over orders of mag- 
nitude in distance for both cases). The most active star formation 
is naturally associated with the deepest part of the potential well, 
so that these results can be interpreted as confirming the intuitive 
picture of Figure [T] 

In both the QSO and GRB case there is a tail of impact pa- 
rameters extending to 100 kpc distances, arising in high mass ha- 
los (A/vir > IO^^Mq) which are dynamically young and contain 
multiple H I clouds hosting active star forming regions. Thus care 
is needed in interpreting the highest 6- values; even though in each 
case the brightest star forming regions are associated with the po- 
tential minimum, observers could expect to observe significant star 
formation closer to the GRB sky coordinates than these extreme 
values would at first suggest. In addition, small halos have such 
slow star formation that they may be invisible for practical pur- 
poses. 

Furthermore, some caution should be exercised in interpreting 
&- values smaller than a kiloparsec, since these lie below the gravita- 
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Figure 3. The distribution of impact parameters b for GRB-DLAs (solid 
thick line) and QSO-DLAs (dotted thick line). The GRB-DLA population 
typically originate very close to the centre of the dark matter halo, whereas 
the QSO-DLAs tend to pass several kpc from the centre, although in both 
cases there is a long tail to large impact parameters. This tail is associ- 
ated with massive (Afvir > 10^^ Mq) halos with recent merging and 
hence multiple gas/star formation centres. The longitudinal distances. A, 
over which the absorbers are spread (see text for exact definition) are also 
shown as thinner lines; the GRB-DLAs are seen to be more compact along 
the line-of-sight. All units are physical. 



Figure 4. The differential column density distribution of GRBs as observed 
(points with error bars; see text for references) and simulated (solid line); 
these are in agreement except for systems with Afjji < 10^^ cm~^. The 
dotted line shows the results obtained if stellar ionizing sources are ignored, 
leading to a slight overestimate of very high column density systems. The 
dashed line shows the results expected from a naive argument based on 
the observed QSO-DLA distribution combined with the Kennicutt star for- 
mation relation; this underestimates the incidence of high column density 
systems, which may be a useful result for understanding the high redshift 
ISM (see text for details). 



tional softening scale of the coarsest simulation included (Table[TJ. 
Reassuringly, however, we found that the only effect of removing 
the low resolution C50 box from our sample was to truncate the 
fe- value distribution above approximately 20 kpc, simply reflecting 
the removal of all high mass halos from our sample: no change in 
the low-6 tail was produced. 

To compare the gas distribution along the sightline of in- 
trinsic and intervening absorbers, we define A to be the distance 
over which the central half of the Hi column is assembled (i.e. 
the physical distance between lower and upper quartile positions). 
The distribution of this quantity is overplotted on Figure |3] as thin 
lines (solid for GRBs, dashed for QSOs). A values suffer equally 
from the resolution caveat mentioned for b above; however we can 
nonetheless conclude that GRB-DLAs are typically assembled over 
a longitudinal distance somewhat less than a kiloparsec whereas 
QSO-DLAs are built up over larger distances (~ 2 kpc). This indi- 
cates that the H I columns in our typical GRB-DLAs are not assem- 
bled from the immediate environment of the GRB itself, but rather 
from gas spread through kiloparsec-sized surrounding regions. 

Although the distributions of b and A are seen to be simi- 
lar to each other, we found that for individual absorbers these two 
quantities are in fact almost uncorrelated; i.e. there is no significant 
relation between the longitudinal dimensions of an absorber and its 
proximity to the central regions of the host protogalaxy. 

3.2 Column density distribution 

For systems at 2; > 1.6 the Lya absorption line is redshifted away 
from otherwise inaccessible UV wavelengths. For such systems, 
it is possible to determine observationally the column density of 



neutral hydrogen. We compiled a list of A^'ni-values for 38 GRB- 
DLAs in the redshift interval 2 < 2 < 4 by combining data from 
[Jakobsson et al. (2006) , |Chen et al.| ( (2007) l and |Fynbo et"ari ( |2009l l. 
We caution that, unlike for the QSO-DLAs for which homogeneous 
catalogues are available (e.g. |Prochaska et al.|2005| l, our GRB-DLA 
statistics are formed from a heterogeneous sample. It is possible 
that extinction of UV/optical afterglows leads to some preference, 
in the observed data, to select lower column density absorbers (with 
correspondingly smaller dust columns; e.g. |Fynbo et al.| 2009); at 
this stage, however, the observational significance of this effect is 
not transparent - further discussion is given in Section [4!4l 

Figure |4] compares the observed column density distribution 
(points with error bars) with our simulated GRB-DLA column den- 
sities (thick solid line). The form of the differential distribution 
is well matched by our simulations for A'^hi > 10^^ cm~^. It is 
immediately clear that the GRB A'^hi distribution is very different 
from the intervening case: for instance, the TVhi distribution di- 
verges at low A^Hi for intervening DLAs but approaches zero for 
GRB-DLAs; further the median GRB-DLA has a column density 
of iVHi ^ 10^^'^ cm"^ whereas only 1% of QSO-DLAs are so 
extreme in a catalogue derived from SDSS DR5 (Prochaska et al.| 
2005 We will discuss these differences and their physical origins 
in Section[3H 

Overall, the simulated column density distribution is in good 
agreement with the observed result. The median A^hi value (TVhi — 

^ The DR5 catalogue is available from 

www.ucolick.org/~xavier/SDSSDLA/. Recently the DR7 
catalogue was analysed with an independent pipeline jNoterdaeme et al"] 
2009 ); the differences are, however, minor for our present purposes. 
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lO'^^'^ cm~^) is reproduced quantitatively. The dotted line in Fig- 
ure |4] shows the output from our simulations without account- 
ing for the effects of local sources (Section |2.1^ ; in this case a 
factor two overabundance of very high column density sightlines 
(A^Hi > W^^'^ cm~^) is predicted. Our simple approach to ac- 
counting for the stellar UV emission has therefore mitigated this 
problem. However, our simulations produce a statistically signifi- 
cant shortfall of systems with A^'hi < 10^^ cm~'^ compared with 
the observed population; we will now discuss an estimate of the 
escape fraction from our simulations, which will highlight this dif- 
ficulty. 



3.3 Escape fraction 

Since we have constructed a sample of neutral hydrogen column 
densities emerging from star formation regions, we can calculate 
the fraction of radiation at the Lyman limit which escapes relative 
to the total produced: 



dNmp{Nm)e 



- CTL Niii 



(4) 



where (Tl — 6.28x10"^* cm^ is the photoionization cross-section 
of a hydrogen atom at the Lyman limit. In effect /esc obeys almost 
identically /esc = p{< 1-6 x lO^'^ cm~^), since the exponential in 
the integral acts as a step function in log A^hi space; however we 
evaluate the full integral since it is not a difficult computation. We 
ignore the effect of dust; although its rate of UV photon absorption 
can be comparable to that of the hydrogen in certain physical situa- 
tions, |Gnedin et al.|p008) argue plausibly that this has little impact 
on the escape fraction since most radiation leaks out along a small 
fraction of extremely low column density (essentially unobscured) 
sightlines. 

Our calculated escape fraction, /esc = 1.0%, is somewhat 
smaller than recent estimates have suggested. In particular, it falls 
short of the 2% ± 0.3% derived by Chen et al. ( 2007 ) from reason- 
able extrapolation of the observed GRB-DLA statistics (see also 
|Fynbo et al.|2009t ; this is a reflection of our previously described 
underestimate of the abundance of absorbers in the observational 
bin 10^'' < Nm/cm-'^ < 10^^ More worryingly, our derived 
fraction is very sensitive to assumptions we make about launching 
the GRB sightlines; for instance, adding a 200 pc gaussian random 
offset to our stellar particle loci raises the fraction to /esc = 2%. 
Such an approach might be justifiable if runaway massive stars are 
preferentially involved in GRB events (e.g. |Hammer et al.|2006^ . 

While some adjustments can increase the reported escape frac- 
tion, others cause it to decrease. For instance when we restricted 
sightlines to originate from younger stellar populations (< 10^ yr) 
than our fiducial technique imposes (Section \2.2\ we derived an 
escape fraction of 0.7%. We also implemented a model in which 
the GRB launch rate is made proportional to the instantaneous 
star formation rate as determined by the gas properties (i.e. tak- 
ing PGRB oc p^'^ in regions where star formation is permitted, as 
opposed to probing new star particles formed by the simulation, 
Section [23| . In this case the escape fraction dropped to just 0.1%. 
These changes appear to be because our stellar feedback algorithm 
rapidly decreases the density in regions around young stars in the 
live simulation, an effect underestimated in the fully instantaneous 
model - but perhaps overestimated in our standard approach which 
allows for excessive delays between formation and GRB events. 

None of our other results - including the column density dis- 
tribution for A^'hi > 10^" cm'^ - are strongly affected by these 



changes, so while the escape fraction seems to highlight shortcom- 
ings in our approximate schemes, we are confident that the majority 
of our results are robust. Further discussion of this point is offered 
in Section|4] 



3.4 Interpretation of the iiigli Nm distribution 

As described above, both observed and simulated GRB-DLAs have 
an entirely different A'^hi -distribution from that of QSO-DLAs. To 
account for this we should consider the relative importance of 
the two major distinctions between these populations (halo mass 
weighting and internal sightline positioning; Section [3T| . We find 
that the difference in the column density distributions arises chiefly 
from the GRB's propensity to probe inner regions of the halo, rather 
than any difference in the column density distribution as a function 
of halo mass. 

We verified this by (unphysically) weighting our QSO-DLA 
sightlines using the GRB-DLA halo weights (wh', equation [3]l - 
the result was an increase in the median A'^hi of just 0.1 dex, 
from 10^°-^ to lO^" ** cm"^. Conversely, by weighting our GRB- 
DLA sightlines using the QSO-DLA halo weights (based on the 
Hi cross-section) we found that the median GRB-DLA column 
density decreased by just 0.1 dex. This result was foreshadowed 
by figure 8 of P08 which shows that the mass of the hosting halo 
has relatively little effect on the A'^hi value of a typical sightline. 

The combination of an observationally-derived fit to the QSO- 
DLA column density distribution (e.g. [Prochaska et al.|2005^ and 
the Schmidt-Kennicutt relation betwee n gas surface den sity and 



Kennicutt 



1998b I actually 



star formation rate (E* oc N^i,n ~ 1.5; 
predicts a GRB-DLA column density distribution. This prediction 
is independent of the details of the underlying halo populations, 
relying only on the QSO absorbers to trace the unbiased neutral 
hydrogen distribution. If the column density distribution follows 
a Schechter function ( |Pei & FaII| 1995] [Prochaska et al.|2005i ) with 
asymptotic slope a and exponential scale A^^, one has the following 
result for the distribution of GRB A'hi values: 



dP 



exp 



^2iVHi 



(5) 



where the 1 in the exponent arises from the change of variables 
A'hi ^ logio -^hi and the 2 in the exponential arises because, by 
symmetry, a sightline to the GRB intersects on average half of the 
total gas column density of the host galaxy. 

We exhibit this function as a dashed l ine in Figurep] using the 
peak posterior dust-correctecj^parameters i Pontzen & Pettini 2009 1 
determined for the previously mentioned SDSS DR5 sample in the 
redshift interval 2.2 < z < 4: a = -1.79, A^.^ = lO^^ '^^ cm"^ . 
Qualitatively, this approach predicts the correct form for the GRB- 
DLA A'^Hi distribution, with a peak and long tail to low column 
densities. However, quantitative details at high A'hi are in some 
tension with the observed GRB-DLA results. (One is less interested 
in pursuing exact agreement below A'^hi < 10^"'^ cm~^ since such 
systems are not included in the fit to the quasar data.) 

P08 showed that our simulated QSO-DLA column density dis- 
tribution, while in good overall agreement with the observed distri- 
bution, over-predicts the occurrence of high A''hi quasar systems. 
Yet the present considerations show that the observed QSO-DLA 



The dust correction is small and does not in fact have a qualitative impact 
on the following considerations, but we include it to remove suspicion that 
dust biasing of QSO samples could account for any discrepancies. 
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distribution Mnder-predicts thie occurrence of high A^hi GRB sys- 
tems - while the simulations are in better agreement with obser- 
vations! This three-way tension between the observed QSO-DLAs, 
observed GRB-DLAs and simulations may be a fruitful line of fu- 
ture exploration to understand details of the high redshift ISM. For 
instance, high-A^Hi quasar absorbers are plausibly limited by H2 
formation l |Schaye|2001^ . The relative over-abundance of high-TVei 
absorbers in GRB afterglows might then suggest dissociation of 
molecular clouds in the vicinity of UV-bright young stellar popu- 
lations responsible for the GRB events. Given the crudeness of our 
gas phase modelling, we are not in a position to make any certain 
conclusions at this stage. 

We should take a moment to consider whether it is appropri- 
ate to apply the Kennicutt relation in the sense above. In a direct 
test of the the relation's applicability to quasar DLA populations, 
|Wolfe & Chen| ( (2006[ > used the A^hi distribution to predict a back- 
ground level of faint diffuse emission at rest-frame wavelengths 
A ~ 1500 A. After blind searches in deep imaging, they con- 
cluded that the Kennicutt relation was over-estimating star forma- 
tion in DLAs. However the authors also noted that this result may 
be negated if a substantial fraction of DLAs are associated with 
nearby bright compact sources, which were excluded from their 
search algorithms. Since this association indeed appears to hold in 
our simulations (P08), the relation may tentatively be applied. We 
note that the naive Kennicutt weighting applied to the simulated 
QSO-DLAs does produce the approximate distribution of simu- 
lated GRB-DLAs, so this method cannot constitute an entirely spu- 
rious comparison. |Nagamine et al.| ( |2008) also placed confidence 
in such a weighting, using it to adapt their simulated QSO-DLA 
statistics to discuss the expected redshift evolution of the GRB- 
DLA A'hi distribution. 
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Figure 5. The neutral gas metallicity of our simulated GRB-DLA sight- 
lines (solid line) compared with the same diagnostic of QSO-DLA sight- 
lines (dashed line). GRB-DLAs have signiiicantly higher metallicities as 
they are associated with star forming regions of more massive halos than 
QSO-DLAs. This picture is broadly supported by the handful of measure- 
ments of (and lower limits on) GRB neutral gas metallicities at 2 > 2, the 
range of which is shown by the horizontal bar. The shaded region shows 
the range of the median if we assume systems with line saturation have 
up to three times the metallicity of their measured lower limits (see text 
for details). The dotted line shows the simulated results after imposing a 
restriction that GRBs can only arise in stellar progenitors with metallicity 
Z < Zq/3. 



3.5 Metallicity Distribution 

Each of our sightlines is assigned a neutral gas metallicity accord- 
ing to equation |2](; using our fiducial techniques we can then pre- 
dict the metallicities observed in cosmological GRB-DLA samples. 
Figure |5] shows the resulting distribution as a solid thick line, with 
the median indicated by a short vertical line bisecting the distri- 
bution. Our metallicity distribution peaks at around Zq/3, with a 
sharp decline to higher metallicities (and long tail to lower); as a 
result the median is approximately Zq /6. 

As a rough observational comparison, we used the lists of 
[Prochaska efaLl ( |2007l |2009) and [Ledoux eTaE] l |2009t to com- 
pile eight GRB-DLA metallicity measurements and a further seven 
lower limits at 2; > 2; to our knowledge this is a complete census 
of such constraints. The small numbers, line saturation effects and 
unknown biase^make it hard, at present, to attempt a quantitative 
comparison of the observed and simulated metallicity distributions. 
We merely indicate in Figure|5]the range of the observed metallici- 
ties (by a horizontal bar), with a vertical line indicating the median 
value. The shaded region shows the range of the median if we in- 
clude assumed metallicity measurements of up to three times the 
value of the lower limits. 

Given the uncertainties our simulations seem to provide an 
acceptable match. Although the observational constraints are not 



tight, there is no doubt that, for instance, our QSO-DLA metallic- 
ity distribution (shown as a dashed line with its median similarly 
indicated) would be deemed inconsistent with the dat^ 

Unlike for the column density distribution, we attribute much 
of this metallicity shift to the 1 dcx higher halo masses probed by 
the typical GRB-DLA compared with a QSO-DLA (Section [STTJ, 
in agreement with [Fynbo et al.| (j2008j. Mean neutral gas metal- 
licity in our halos range from ~ 10"'^ Zq for M^ir ~ W^Mq 
to ~ IO^^Zq for i\fvir = W^^Mq: a 0.7dex rise is expected 
from this effect. The remainder of the effect is accounted for by 
the marginally higher metallicities of the central regions which the 
GRB sightlines probe. This is a small ^ O.ldex correction, de- 
spite large internal abundance gradients discussed below, because 
the sightline averaging smears out much of the metal density vari- 
ations. 



3.6 Effect of GRB Metallicity Ceiling 

So far we have assumed that GRB events simply trace, in an unbi- 
ased fashion, star formation activity in the cosmos. However, GRBs 
are thought to arise from rapidly rotating progenitor stars. Since an- 
gular momentum is concentrated in outer layers of high mass stars, 
it is easily lost through winds. Mass loss rates are thought to scale 



^ If anything, the observed sample is likely biased in favour of low metal- 
licity GRB-DLAs which have smaller dust columns and therefore brighter 
optical afterglows; see e.g. |Fynbo et al.|j2008| . 



^ Our new simulations have a QSO-DLA metallicity distribution with a 
median [Z/Zq] ~ —1.6, approximately 0.2dox lower than the median in 
P08, as described in Section|2] 



© 2009 RAS, MNRAS 000, 000-000 



Simulating GRB afterglow absorbers 9 




-2 2 

x/kpc (physical) 



-2.0 



-1.5 



-1.0 



-0.5 



0.0 



Figure 6. Map of the metallicity of the gas in a thin slice through the centre 
of the most massive halo in our MG box (Afvir ~ 3 X lO^^M©), with the 
volume density of H I overplotted as contours (the outermost white contour 
corresponds to rim = 10~^ cm~^ and inner, black, contours indicate 
0.1, 1, 10 and 10^ cm""'). Little of the gas by volume exceeds 1/3 solar 
(red region, the cut-off in our truncated GRB model), and star formation in 
nearby regions can still give rise to GRBs with lines-of-sight intersecting 
the high Z region. 



almost linearly with metallicity, so GRBs could well arise preferen- 
tially from low metallicity progenitor stellar populations ^Woosley] 
|& Bloom|2006| and references therein). There is also some obser- 
vational evidence, particularly at low redshifts, that GRB-hosting 
populations are indeed biased in this sense (e.g. [Fruchter et al 
2006, 



Modjaz et al] 2 008, [Calura et al.| 2009; although see |Wolf 



|& Podsiadlowski| 2007, |Kelly et al.| 2008, |Savaglio et"aL] 2009). 

We tested the effect of such a picture on our results by pro- 
ducing a new set of sightlines, allowing GRBs to be launched 
only from stellar populations with Z < Zq / 3 (as suggested by 
[Woosley & Heger||2006[ > and re-calculating the halo weights ac- 
cording to w'f^ — Whf{Z < ZqI?>) where Wh is defined by and 
f{Z< Zq /3) is the mass fraction of young stars (according to our 
criterion in Section|2l2]l with metallicity below the adopted limit. 

We found that the differences in column density and impact 
parameter distributions are entirely negligible; the only significant 
effect was seen in the sightline neutral gas metallicities and the stel- 
lar population metallicities (the latter quantity is discussed below 
in Section [3J) l. The updated neutral gas metallicity distribution is 
shown as a dotted line in Figure[5]and shows a reduction in median 
metallicity by 0.2 dex to [Z/Zq] = —0.9. This is a relatively small 
shift, largely reflecting that even our most massive halos at z = 3 
have mean cold gas metallicities of only Zq /2. 

Furthermore sightlines of up to solar metallicity remain at- 
tainable with the cut because the neutral gas within each halo has 
metallicity dispersions of up to an order of magnitude; thus sight- 
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Figure 7. Properties of the stellar population within the host galaxies of 
z = 3 GRB (thick solid line, dotted line for < Zq / 3 progenitors) and 
QSO-DLAs (dashed line). From top to bottom we plot the probability den- 
sities distributions for star formation rate, total stellar mass and (mass- 
weighted) stellar population metallicity respectively. A second scale on the 
uppermost panel converts the star formation rate to an ultraviolet AB mag- 
nitude as described in the text, ignoring dust obscuration. The properties of 
the host galaxies of the two classes of absorber differ in accordance with 
their different virial masses, with the GRBs favouring somewhat higher star 
formation rates, total stellar masses and stellar metallicities. This applies 
even for the cut metallicity GRB population, although the high-Z peak is 
slightly suppressed in this case. 



lines which start in a low metallicity environment may nonetheless 
encounter a nearby dense knot of solar metallicity gas which dom- 
inates the line-of-sight measurement. This is illustrated by Figure 
|6] which shows the metallicity of the gas in a thin slice through 
the centre of the largest halo in the MG box. Red and white regions 
have metallicity such that they are excluded from hosting GRB pro- 
genitors. Overplotted contours show the regions in which the neu- 
tral gas volume density exceeds nni = 10~^, 0.1, ■ ■ ■, 10^ cm""^ 



(respectively from outer to innermost contour). Dense neutral re- 
gions (> 0.1 cm^'') have mean metallicity lO"" * Zq with a dis- 
persion of 0.8 dex. Qualitatively GRBs are to be found in the dens- 
est regions which are not forbidden by the metallicity cut (i.e. not 
within the red area); the prime location remains very close to the 
centre so that high metallicity gas will present a large cross-section 
to the sightline. 



3.7 GRB host halo properties 

Our simulation naturally associates each QSO- or GRB- sightline 
with a host halo, the properties of which can readily be computed. 
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In this section we briefly describe the distribution of star formation 
rates, total stellar masses and stellar metallicities in the populations 
within each halo. 

Figure|7] shows the distribution of each of these properties (in 
three separate panels) with the GRB-DLAs and QSO-DLAs repre- 
sented by solid and dashed lines respectively. The dotted lines show 
the population of GRB-DLAs if a progenitor metallicity criterion 



Z < Zq/?, is introduced (see Section 3.6 1. The GRB-DLA hosts 
generally show somewhat higher star formation rates, total stellar 
masses and metallicities than the intervening population: mean star 
formation rate^for the GRB and QSO DLAs are 10 M© yr"^ and 
lM0yr~^, a difference of Idex. However this should be com- 
pared with a wide dispersion of around 4dex in both cases. We 
have added to the top panel a scale showing corresponding unob- 
scured UV magnitudes using the luminosity conversion ratio from 
|Kennicutt| ( [ 1 998a^ , corrected by a factor 1 .6 in order to consistently 
use the Kroupa IMF assumed by our simulations. 

Mean stellar masses for the QSO-DLA and GRB-DLA host 
populations differ by around an order of magnitude (with respec- 
tive values 10*'^ Mq and lO^'-^M©), but exhibit wide disper- 
sions around these values. For the stellar metallicity (distinct from 
the neutral gas metallicities described in Section |3.5^ the mean 
values are 10"° '' ^0, 10'° '-' for unbiased and metallicity- 
biased GRBs respectively compared against 10^^'* Zq for the 
QSO DLAs. 

The agreement of these results with observational constraints 
have been discussed elsewhere for QSOs (P08). For GRBs, no ho- 
mogeneous high redshift catalogues of host galaxies have been as- 



sembled. However the four 2: J> 1 objects studied by Christensen 
|et al.|p004| have star formation rates well within the range sug- 
gested by our simulations; see also [Jakobsson et al.| ( |2005^ . At 
2 > 2, the eight GRB-DLAs with identified host galaxies listed by 
|Fruchter et al.| ( |2006) have a range of luminosities — 16 ^ Mv 1? 
—22 in the observer-frame V band; broadly, given that our brightest 
objects should be substantially dust obscured in the rest-frame UV, 
our results are consistent with these observations. 

All of the distributions can be interpreted by inspecting the 
relation between virial mass and the property under consideration 
(see, for example, figure 13 of P08 for the star formation rate). As 
well as satisfactorily accounting for the approximate means in each 
case, specific features in the distribution function can be explained 
in this way. In particular Z(Mvir) rises steeply from 10"'^ Zq at 



10^ Mq to around half solar at lO^M©, beyond which the 



metallicity remains constant. Physically this is because massive ha- 
los harbour older stellar populations as well as the young enriched 
stars from recent formation activity. This gives rise to bimodality 
in the metallicity of the QSO-DLA stellar populations, with peaks 



at Z 



10' 



Zq (corresponding to the metallicity of the peak 



QSO-DLA mass, 10^" Mq) and Z ~ lO""'^ Zq (corresponding 
to the plateau metallicity: all halos with M > lO" Mq fall into 
this bin). Bimodality is not seen in the neutral gas metallicities be- 
cause these measurements are averaged only along the line of sight, 
not the entire halo, and local metallicities can reach values of up to 
solar: see Figure [6] 



* Calculated as described in P08 section 4.3. The QSO distribution shown 
in the middle panel of Figure|7]is equivalent to figure 14 of P08; however the 
latter plot shows dp/dAfj,, not dp/d logjg M^, as its label in the published 
version erroneously claims - this was a typographical error 



4 CONCLUSIONS AND DISCUSSION 
4.1 Overview 

We have investigated the properties of high redshift (z ~ 3) GRB- 
DLA absorbers in a series of galaxy formation simulations, on the 
assumption that long duration GRBs are associated with the death 
of massive, short-lived stars. Following our success in matching 
properties of traditional DLAs along randomly chosen sightlines to 
quasars (P08), we used our results to interpret observational differ- 
ences between the populations, finding two primary distinctions. 

Firstly, because the star formation rate scales more steeply 
with virial mass (Af* ~ AfJif,P08fi gure 13) than does the neutral 
hydrogen cross-section (approximately (Tdla ~ M^ir, although a 
single power-law is an oversimplification; P08 figure 4), the halos 
contributing most to the GRB-DLA incidence rate are considerably 
more massive than those contributing to the QSO-DLA incidence 
rate (Section [3T| l. The dominant contribution for GRB-DLAs arises 
in halos with 10^° < Afvir /Mq < lO"; for QSO-DLAs this 
range is an order of magnitude lower, 10^ < MvIt/Mq < 10^^. 
While this argument obviously relies on our simulated scalings 
being accurate, there are many reasons for being cautiously opti- 
mistic: in particular, the simultaneous matching of the QSO-DLA 
metallicity distribution (P08) and the galaxy mass-metallicity rela- 
tion ^Brooks et al.|2007 1 is hard to reproduce unless the underlying 
relations are realistic. |Fynbo et al.H2008[ > have previously reported 
that GRB metallicities can be accounted for if a Holmberg-like re- 
lation, ctdla oc Af^* with t ~ 0.4, holds at high redshift; here 
we have shown that the metallicities are reproduced and that just 
such a scaling arises naturally from cosmological simulations (al- 
beit with t ~ 0.3). The difference in mass ranges has some effect 
on the predicted properties of the hosting halos, most notably lift- 
ing the mean stellar metallicity from lO"^ * Zq for QSO-DLAs to 
10"*^'^ Zq for GRB-DLAs. However, these effects are relatively 
small compared to the large spread of hosting halo properties in 
both populations, reflecting that incidence rates for both classes are 
a shallow function of virial mass. 

[Courty et al.| ( |2004[ l took a reciprocal approach, using ob- 
served galaxy properties to determine the likely underlying popula- 
tion of host objects. Their results showed that GRB hosts at 2 = 3 
have typical masses A/vir < IO^^Mq which does not agree with 
our analysis. This may be because [Courty et aL|p004)) adopt, on 
an observational basis, the star formation efficiency (rate of star 
formation per unit stellar mass) as an indicator of GRB probability. 
Our simulations do not suggest a strong link between star formation 
efficiency defined in this way and the GRB event rate; if extended 
observational analysis continues to point to such a link, further in- 
vestigation will be required. 

The second distinguishing effect we find is that the GRB- 
DLAs arise in gas closer to the centres of the hosting objects than 
a conventional DLA (as seen observationally, e.g. [Savaglio et al.| 
|2003||Vreeswijk et al.|2004[|Fruchter et al.|2006^ . While resolution 
artifacts cloud an exact interpretation of some details, we are con- 
fident in concluding that approximately 50% of GRB-DLAs arise 
less than 1 kpc from the centre of their parent halo, whereas this fig- 
ure for QSO-DLAs is just 10% (and, accordingly, the QSO-DLAs 
have a higher median impact parameter of ~ 4 kpc). Similarly, 
GRB-DLAs tend to be more concentrated along the line of sight; 
we defined the measure A to be the distance over which half the 
total hydrogen column density is accumulated and found that this 
took a median value of 1 kpc for GRB-DLAs and 2 kpc for QSO- 
DLAs. This indicates that the neutral hydrogen seen in GRB af- 
terglow spectra is not to be found in the immediate circumstellar 
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environment, but is spread over an extended region - although a 
detailed analysis would require better resolution of the ISM than 
our present simulations are able to provide. 

4.2 Observational diagnostics 

We investigated the neutral hydrogen column density (Sec- 
tions |3.2| |3.4^ and the metallicity (Section |3.5^ distributions of 
our simulated absorbers, finding good observational agreement in 
both cases except for a deficiency in low column density (A*'hi < 
10^^ cm^^) absorbers. The simulated neutral hydrogen column 
densities, as for the observed GRB absorbers, are much higher 
than those seen in intervening DLAs. We conclude that this is 
largely due to the lowered impact parameters, with the change in 
the masses of halos having a secondary effect (for instance, figure 
8 of P08 shows that the relation between column density and halo 
mass is relatively weak). 

The metallicities are higher, by almost a factor of ten, than 
those seen in QSO-DLAs - even if GRB-DLAs can only arise from 
stellar populations with Z < Zq /3 (a limit motivated by the col- 
lapsar picture, e.g. |Hirschi et aLpOO?) |Woosley & Heger||2006| l. 
This suggests that the observed population of high redshift GRB- 
DLAs and GRB hosts do not pose any severe problems for the typ- 
ical coUapsar GRB paradigms, counter to the conclusions of |Wolf| 
|& P odsiadlowski' (2007^ and' Lapi et"S^ ( [2008| >, but in agreement 
with Nuza et al. ( 2007 ). The ineffectiveness of a metallicity ceiling 
on the GRB-DLA progenitors to limit the metallicity of the overall 
sightline population arises because, within each protogalaxy, neu- 
tral gas metallicity dispersions of up to an order of magnitude are 
possible (greatly exceeding the 0.2 dex allowed by |Wolf & Pod-| 
|siadlows ki]2007). In our picture GRBs now originate in dense re- 
gions with lower metallicities, but they are typically still close to 
gas with significantly higher metallicity (Figure [6| which conse- 
quently presents a large covering fraction. Furthermore high mass, 
metal rich halos continue to contribute significantly to the GRB 
rate: even for our most massive galaxy at 2; = 3 (with a mean neu- 
tral gas metallicty of ~ Zq/2), 30% of the star-forming gas mass 
remains at low metallicities Z < Zq/3. We note that, since our 
simulations include a turbulent metal diffusion term, this disper- 
sion cannot be attributed solely to numerical effects (Section|2j. 

The leaves open the question of whether the observational sit- 
uation at lower redshifts (z < 1.5), where host galaxies may be 

ffruchter et al.| 2006; 



less massive than typical hosts of supernovae 



although see |Kelly et al.| 2008, |Savaglio et al. 2009), can be quan- 
titatively reconciled with less pronounced differences between su- 
pernovae and GRB hosts at high redshift. A natural explanation is 
simply that at earlier times mean metallicities are lower and high 
metallicity gas is typically less well mixed W uza et al.|2007[|Ko-| 
[cevski et al.|2009} . In principle we could probe this claim in our 
simulations, but such work is beyond the scope of the present in- 
vestigation. 



4.3 UV Escape fraction 

One of the shortcomings of our work is the unreliability of our pre- 
dicted 1% escape fraction of UV photons from gamma ray burst 
regions. This is somewhat lower than recent simulations of escape 
fractions (~ 2% from results by|Gnedin et dl]2008; between 2% 
and 10% according to |Razoumov & Sonuner-Larsen 2006; or, for 
halos with masses M < lO"M0,|Wise & Cen|2009 reported es- 



cape fractions between 25% and 80% in very high resolution Eule- 
rian simulations). Furthermore it falls short of constraints on escape 



fractions from comparing rest-frame UV emission on either side of 
the Lyman limit break (~ 2%, although this value is sensitive to as- 
sumptions about the dust extinction; |Shapley et al.|2006) and most 
directly underpredicts the escape fraction implied by current GRB 
observations (~ 2%; jChen et al.|2007 also visible directly in Fig- 
ure[4]l. However, most importantly, it is not robust to changes in our 
sightline construction algorithms (Section [3.2^ , unlike our other re- 
ported results. 

The fragility reflects enormous uncertainties in this regime. In 
addition to the inherent complexities in accurately simulating ab- 
sorption systems, for GRB-DLAs one can no longer justify the as- 
sumption that local ionizing UV sources are unimportant compared 
with cosmological sources (Section [2.1^ . We have taken a simple 
pragmatic approach to this problem by noting that the ionizing ra- 
diation of local sources is largely absorbed by H II regions which 
we paint into our output. (More accurate radiative transfer would 
probably lend confidence to our results, but we caution that with- 
out properly resolved Hll regions the physical meaning of even 
the most precise transfer algorithms is still obscure.) Furthermore, 
there is a sensitivity to the energy injection from supernovae, as 
demonstrated by the drop in escape fraction to 0.1% when trigger- 
ing GRBs based on the star formation rate scaling of our gas par- 
ticles (rather than existing young stellar particles which will have 
injected feedback energy into the nearby ISM). Finally, if runaway 
stars are important contributors to the GRB incidence rate l |Han>| 
|mer et al.|2006| l one could expect a significant increase in the num- 
ber of low column density absorbers. 



4.4 Model Uncertainties 

There are further complications specific to GRB-DLAs which we 
have not specifically addressed. One is the evidence, from fine- 
structure transitions, that the neutral phase metals detected in the 
systems are found at distances ranging from hundreds of par- 
sees (Prochaska et al. 2006b to as much as a kiloparsec or more 
( [Vreeswijk et al..2007 , D'Elia et al.|2 009). It is possible that the 
GRB afterglow itself ionizes a significant amount of material along 
the line of sight. Given the known afterglow curves, the ioniza- 
tion of material by the GRB itself has been advanced as a pos- 
sible explanation (e.g. [Vreeswijk et al.| 2007), although whether 
bubbles of more than a few parsecs radius can be ionized is un- 
clear (e.g. [Prochaska et"ar] 2008). The afterglow-ionizing scenario 
is certainly not firmly established: perhaps the hot progenitor stellar 
populations have already carved out large ionized bubbles within 
the ISM, or the densest knots of gas which are metal rich are sim- 
ply found at some distance from the GRB site (as we have men- 
tioned above). Unfortunately we do not have resolution to probe 
the required scales reliably, but we did produce comparison statis- 
tics in which we ionized all material within a 1 kpc sphere of a 
GRB site before measuring its sightline properties. This led to a 
decrease in the numbers of very high column density absorbers 
{Nhi > 10^^ cm"^), reducing the median to A^m = lO^'^-^cm^^ 
While the agreement of our simulations with observations was 
therefore impaired, our results were not qualitatively altered. Simu- 
lations with better resolution of the ISM will be necessary to prop- 
erly constrain the effects of GRBs on their immediate environment. 

Similarly we have been unable to make a quantitative assess- 
ment of the possible biasing effects of dust. It is possible that dust 
obscuration of optical afterglows significantly biases against ob- 
serving high metallicity, high column density absorbers. Indeed the 
phenomena of 'dark bursts', GRBs for which the optical afterglow 
is not seen or is dimmed by several magnitudes, could be explained 
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if these bursts are simply obscured by large columns of dust ( |Perley| 
|et al.|2009p . However, it appears to be inappropriate to extrapolate 
the dust columns predicted for quasar absorbers (by models such 
as those in |Pontzen & Pettim|2009) to our simulated GRB results: 
such an approach would significantly over-predict the reddening 
seen in real GRB-DLA samples (for a compilation see |Fynbo et al.| 
\2QQ9) and would suggest, for example, th at GRB 080607 (with 

2009) 



Prochaska et al. 



A'^Hi = 10 cm"^, Z ^ Zq according to|] 
should have been completely dark. Either the GRB-DLAs arise in 
regions in which simple scalings of dust columns fail, or perhaps 
the intense outgoing radiation is able to rapidly destroy dust even 
beyond the hydrogen ionization front. Whatever the explanation, 
we have been unable to make full sense of this issue in our present 
exploration and leave its resolution to future work. 

Despite these difficulties, we have affirmed in the present work 
that GRB-DLAs offer a new insight into the inner regions of the 
protogalactic population at high redshifts. Their basic statistics can 
be matched by simple techniques applied to existing cosmological 
simulations. There can be little doubt that, as observational datasets 
increase in size and sophistication, GRBs and their associated ab- 
sorbers will make a significant contribution to our understanding of 
the structure and formation of galaxies. 
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